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ABSTRACT 

We report the discovery of the galaxy cluster C1J1226.9+3332 in the Wide Angle ROSAT Pointed 
Survey (WARPS). At z = 0.888 and L x = 1.1 x 10 45 h" 2 erg s" 1 (0.5 - 2.0 keV) C1J1226.9+3332 is the 
most distant X-ray luminous cluster currently known. The mere existence of this system represents a 
huge problem for Qq = 1 world models. 

At the modest (off-axis) resolution of the ROSAT PSPC observation in which the system was detected, 
C1J1226. 9+3332 appears relaxed; an off-axis HRI observation confirms this impression and rules out 
significant contamination from point sources. However, in moderately deep optical images (R and I 
band) the cluster exhibits signs of substructure in its apparent galaxy distribution. A first crude estimate 
of the velocity dispersion of the cluster galaxies based on six redshifts yields a high value of 1650 km s _1 , 
indicative of a very massive cluster and/or the presence of substructure along the line of sight. While 
a more accurate assessment of the dynamical state of this system requires much better data at both 
optical and X-ray wavelengths, the high mass of the cluster has already been unambiguously confirmed 
by a very strong detection of the Sunyaev-Zel'dovich effect in its direction (Joy et al. 2001). 

Using C1J1226.9+3332 and C1J0152. 7-1357 (z = 0.835), the second-most distant X-ray luminous 
cluster currently known and also a WARPS discovery, we obtain a first estimate of the cluster X-ray 
luminosity function at 0.8 < z < 1.4 and L x > 5 x 10 44 h^ 2 erg s -1 (0.5 — 2.0 keV). Using the best 
currently available data, we find the comoving space density of very distant, massive clusters to be in 
excellent agreement with the value measured locally (z < 0.3), and conclude that negative evolution is 
not required at these luminosities out to z ~ 1. Our findings are in conflict with earlier claims of highly 
significant (> 3er) negative evolution already at 0.3 < z < 0.6 based on the cluster samples of the EMSS 
and the CfA 160 degree survey. Our results agree, however, with the lack of significant evolution of very 
X-ray luminous clusters out to z ~ 0.4 reported by the MACS team. Our findings are also consistent 
with the abundance of very X-ray luminous clusters at z ~ 0.8 inferred from the EMSS cluster sample, 
provided C1J0152. 7-1357 (which was missed by the EMSS) is added in. 

Subject headings: galaxies: clusters: general — galaxies: clusters: individual (C1J1226. 9+3332, 
C1J0152. 7-1357) — cosmology: observations — X-rays: general 
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1. INTRODUCTION 

Measurements of the abundance of clusters of galaxies 
as a function of redshift allow a number of physical and 
cosmological parameters of structure formation models to 
be constrained (e.g. Oukbir & Blanchard, 1997; Eke et 
al., 1998). The tightest constraints are obtained from ob- 
servations of the most massive and most distant clusters 
which are extremely rare in all models of cluster formation. 
For instance, the predicted space density of galaxy clus- 
ters with intra-cluster gas temperatures of about 8 x 10 s K 
(kT ~ 7 keV) at z ~ 1 is two orders of magnitudes higher 
in an open universe with f^o = 0.3 than in a flat universe 
with n = 1 (Viana & Liddle, 1996). 

Among the many efforts to compile statistically rele- 
vant samples of these rare systems, X-ray flux limited sur- 



veys carry particular appeal. Galaxy clusters are bright 
X-ray sources that can be detected out to high redshift. 
Moreover, X-ray emission from clusters originates from gas 
trapped and heated in deep gravitational potential wells. 
X-ray surveys thus naturally select three-dimensionally 
bound systems and are almost unaffected by projection ef- 
fects. Finally, X-ray flux limited surveys greatly facilitate 
the measurement of comoving space densities since effec- 
tive search volumes are easily computable from the sur- 
veys' selection functions which are usually a simple func- 
tion of X-ray flux and, for distant clusters, almost inde- 
pendent of X-ray source extent. 

While one would ideally want to study clusters that are 
both distant and massive, the rarity of these objects cur- 
rently forces the observer to give priority to one of the two 
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qualifiers. Previous and ongoing cluster surveys have thus 
adopted either of two fundamentally different, but com- 
plementary, approaches. For instance, very deep X-ray 
and optical surveys covering only small areas of sky have 
become increasingly successful at finding poor clusters at 
z }> 1 (Rosati et al., 1999; Lubin et al., 2000) while, at the 
other extreme, relatively shallow X-ray surveys covering 
very large solid angles (Ebeling, Edge & Henry, 2001) are 
in the process of producing large samples of very massive 
systems at lower redshift (z ~ 0.5). 

2. THE WIDE ANGLE ROSAT POINTED SURVEY 

The Wide Angle ROSAT Pointed Survey (WARPS) is 
one of a few surveys straddling the dividing line between 
the two strategies outlined in the previous section. Fol- 
lowing the approach pioneered by the EINSTEIN Ex- 
tended Medium Sensitivity Survey (EMSS, Gioia et al., 
1990; Stocke et al., 1991), WARPS searches for distant 
clusters among a large number of X-ray sources serendip- 
itously detected in pointed ROSAT PSPC observations. 
The WARPS strategy, as well as earlier results, have been 
discussed in previous papers (Scharf et al., 1997; Jones et 
al., 1998; Ebeling et al., 2000; Fairley et al., 2000). 

3. DISCOVERY OF CLJ1226.9+3332 

C1J1226.9+3332 was detected in the WARPS survey as 
an extended X-ray source 14.5 arcminutes off axis in the 
ROSAT PSPC observation of NGC 4395, a nearby, low- 
luminosity Syl galaxy (RP600277N00, on-axis exposure 
time 9036s). The significance of detection exceeds 10 a in 
a 2 arcmin (radius) aperture. The source is also detected 
at more than 5 a significance (within an aperture of 1' 
radius and at an off-axis angle in 14.1') in a ROSAT HRI 
observation of the same target (RH702725N00, on-axis ex- 
posure time 11,353s). 

Figure 1 shows an I band image of the source with adap- 
tively smoothed X-ray contours from the ROSAT PSPC 
and HRI observations overlaid. The emission is well cen- 
tered on an apparent overdensity of faint galaxies. 

3.1. X-ray observations 

The net PSPC count rate directly detected by the VTP 
algorithm (Ebeling & Wiedenmann, 1993; see also Scharf 
et al. 1997 for details of how WARPS employs this de- 
tection algorithm) is (2.29 ± 0.18) x 10~ 2 ct s Ll in PHA 
channels 50 to 200. The exposure time at the location of 
the source (a = 12 h 26 m 58.1 s , 6 = +33° 32' 50.1", J2000) 
is 7.9 ks. Assuming a beta model (Cavaliere & Fusco- 
Femiano, 1976) with (3 = 2/3 and core radius 70 h^ Q kpc, 
derived from the distribution of the immediately detected 
photons (sec Scharf et al. 1997 for details), and taking 
into account the effects of the PSPC off-axis point-spread 
function (PSF) we extrapolate to a total source count rate 
of (2.77 ± 0.21) x 10~ 2 ct s -1 . Using the Galactic value 
of 1.38 x 10 20 cm -2 for the equivalent column density of 
neutral hydrogen in the direction of the source (Dickey & 
Lockman, 1990), a metallicity of 0.3, and a gas temper- 
ature of kT — 11.7 keV (see below) we convert the total 
PSPC count rate into a total X-ray flux in the 0.5-2.0 keV 
band of (3.4 ± 0.3) x 10~ 13 erg s^ 1 cm" 2 . 

Following the same procedure for the HRI data, we find 
an exposure time of 10.2 ks at the location of the source 



(a = 12 h 26 m 58.0 s , £ = +33° 32' 54.1", J2000). This 
source centroid agrees to within 4.2" (less than the typical 
ROSAT astrometry error) with the one determined from 
the PSPC data. The extrapolated total cluster count rate 
is measured to be (1.10 ± 0.18) x 10~ 2 ct s -1 . We use 
the same assumptions as before to convert this count rate 
to a total X-ray flux of (2.9 ± 0.2) x 10~ 13 erg s _1 cm" 2 
(0.5-2.0 keV), a value that is in good agreement with the 
PSPC measurement. Although the HRI image shows no 
obvious point sources within three arcminutes of the source 
(see Fig. 1) we cannot strictly rule out the possibility that 
some or all of the observed X-ray flux originates from a 
concentration of individually faint point sources. 

3.2. Optical observations 

We obtained redshifts of 11 objects close to the appar- 
ent cluster core shown in Fig. 1 using the LRIS spectro- 
graph (Oke et al. 1995) on the Keck-II 10m telescope. 
Both apparent cluster members and possible X-ray con- 
taminants (blue galaxies and stars) were targeted. Com- 
bining the results of a longslit observation carried out 
in January 1999 (300/5000 grism, dispersion 2.6A/pixel, 
slitwidth 1.5 arcsec, wavelength range 5000-10000A) with 
those of multi-object spectroscopy performed in January 
2000 (600/7500 grism, dispersion 1.3A/pixel, slitwidth 1.5 
arcsec, wavelength range typically 6200-8700A) we find no 
obvious contaminants (broad line emitters) but six galax- 
ies with accordant redshifts around a (heliocentric) mean 
of z — 0.8877. The six spectroscopically confirmed cluster 
members are marked in Fig. 1; their positions and redshifts 
are listed in Table 1. The position of the brightest clus- 
ter galaxy (labeled A in Table 1) coincides within 2" (6") 
with the X-ray emission centroid as determined from the 
PSPC (HRI) observations discussed in Section 3T, sup- 
porting our identification of this source as a distant cluster 
of galaxies. 

4. INTRINSIC PROPERTIES OF CLJ1226. 9+3332 

From the six redshifts measured by us so far (see Ta- 
ble 1) we obtain a very crude first estimate of the cluster 
galaxy velocity dispersion of a = 16501^40 km/s. If con- 
firmed in future observations, this very high value suggests 
a very massive cluster and/or substructure along the line 
of sight. 

Assuming a mean cluster redshift of z = 0.8877 we de- 
rive a total X-ray luminosity of L x = (1.06 ± 0.08) x 10 45 
h^ 2 erg s _1 (0.5 — 2.0 keV, go = 0.5 is assumed throughout 
this Letter) from the total cluster flux as determined from 
the PSPC observation, corresponding to an estimated gas 
temperature of 11.7 keV based on the Lx — kT relation 
of White, Jones & Forman (1997). The inferred bolomet- 
ric luminosity of this sytem is 8.0+ 2 ; 2 , x 10 45 h^ 2 erg s _1 
(we assume an uncertainty of ±1 keV in the estimated 
kT value), making CI J1226. 9+3332 more X-ray luminous 
than any other cluster currently known at z > 0.55. The 
cluster's estimated X-ray temperature is consistent with 
its high velocity dispersion (Mushotzky & Scharf, 1997). 

While the X-ray and optical properties of 
C1J1226. 9+3332 described above strongly suggest the 
presence of a highly massive cluster, the limited depth 
of the existing X-ray and optical data does not allow 
us to rule out contamination from a multitude of X-ray 
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faint point sources and/or from projection effects. Un- 
ambiguous evidence of the high mass of this system was, 
however, obtained recently in the form of the detection of 
a very strong Sunyaev-Zeldovich Effect (SZE) decrement 
centered on the X-ray position reported here. The SZE 
measurements were obtained interferometrically at the 
Bcrkclcy-Illinois-Maryland Array (BIMA) at a frequency 
of 28.5 GHz, and are reported by Joy et al. (2001). The 
BIMA observations confirm the presence of a deep gravi- 
tational potential well and lead to an estimate for the total 
gravitational mass of C1J1226. 9+3332 that is similar to, 
and possibly higher than that of the well studied, massive 
cluster MS 1054.4-0321 at z = 0.83 (Hoekstra, Franx & 
Kuijken 2000, and references therein). 

5. THE EFFECTIVE DEPTH OF SERENDIPITOUS X-RAY 

CLUSTER SURVEYS 

All serendipitous X-ray cluster surveys past and present, 
from the EMSS to the ROSAT Deep Cluster Survey 
(RDCS, Rosati et al., 1998), probe deep enough to detect 
very X-ray luminous clusters (Lx > 8 x 10 44 h^ 2 erg s _1 
, 0.5-2.0 keV) out to redshifts z £ 1.5. WARPS, for in- 
stance, would have detected C1J1226. 9+3332 over the full 
geometric solid angle covered (72.0 deg 2 ) out to a maximal 
redshift of z = 1.87. 

Although formally correct, such maximal detection red- 
shifts are misleading. WARPS, like all other serendipitous 
cluster surveys conducted to date (with the notable ex- 
ception of the RDCS), performed its imaging follow-up 
observations exclusively in the optical R and I bands. For 
clusters at z > 1, R and I correspond to B or even U in 
the cluster rest frame and are thus very inefficient bands 
for the detection of distant cluster ellipticals. Indeed, and 
not surprisingly, WARPS 9 failed to detect any cluster at 
z > 1. We thus argue that the limitations of the optical 
imaging observations impose a more stringent constraint 

of ^max ^ 1 • 

Imaging of X-ray selected cluster candidates at near- 
infrared wavelengths (as conducted by the RDCS team) 
allows the identification of clusters beyond z ~ 1 (Rosati 
et al., 1999). However, instrumental limitations still im- 
pose a redshift limit of z ~ 1.4 as, at z > 1.4, the most 
important absorption features (the Ca H and K doublet) 
are redshiftcd to observed wavelengths redward of 9,500 
A where the efficiency of even the most powerful present- 
day spectrograph (LRIS on the Keck 10m telescope) is less 
than 5 per cent. Consequently even the currently deepest 
X-ray cluster survey, the RDCS, has yet to spectroscopi- 
cally confirm a cluster at z > 1.3. 

6. THE EVOLUTION OF THE MOST MASSIVE CLUSTERS 

OF GALAXIES OUT TOZ~l 

We attempt to constrain the evolution of the X-ray 
cluster luminosity function (XLF) at the highest redshifts 
and highest luminosities using the small, but statistically 
complete WARPS sample defined by the selection crite- 
ria z > 0.8 and L x > 5 x 10 44 h^ 2 erg s" 1 (0.5 - 2.0 
keV). The only two WARPS clusters meeting these crite- 
ria, C1J0152.7-1357 (Delia Ceca et al., 2000; Ebeling et al., 
2000) and C1J1226.9+3332, are the most X-ray luminous 

9 as well as all other cluster surveys relying exclusively on imaging 
10 A discussion of the WARPS XLF of less luminous systems at z > 



distant clusters currently known 10 . Both C1J0152. 7-1357 
(which was missed by the EMSS, Ebeling et al., 2000) 
and C1J1226. 9+3332 feature X-ray fluxes well above the 
WARPS flux limit and would have been detected out to 
maximal redshifts of z — 1.57 and z — 1.87, respectively. 
However, for the reasons described in the previous section, 
we use a lower maximal redshift of z — 1 when computing 
the WARPS search volumes. 

The resulting cumulative WARPS X-ray luminosity 
function of very X-ray luminous (Lx > 5 x 10 44 h^ 2 erg 
s" 1 , 0.5 — 2.0 keV) and very distant (z > 0.8) clusters 
is shown in Figure 2. Standard la errors (in the Poisson 
limit, Gehrels 1986) are indicated by the dark, shaded re- 
gion. We also compute an alternative version of the XLF, 
using the very generous assumption of z max = 1.4 (shown 
by the hatched region in Fig. 2). The currently best de- 
termination of the local (z < 0.3) cluster XLF from the 
ROSAT Brightest Cluster Sample (Ebeling et al., 1997) 
is shown by the thick, solid line. We find both WARPS 
XLF estimates to be in very good agreement with the lo- 
cal measurement. Since the search volumes of WARPS and 
EMSS for clusters with L x ~ 9xl0 44 h^ 2 ergs" 1 (0.5-2.0 
keV) at 0.8 < z < 1 are very similar (about 1 x 10^ h^ 
Mpc 3 ) our measurement is also in good agreement with the 
XLF inferred from the two X-ray luminous EMSS clusters 
in this redshift interval, MS1054-0321 (z = 0.829) and 
C1J0152.7-1357 (originally missed by the EMSS). 

Also shown in Fig. 2 is the prediction of the XLF evo- 
lution model derived by Rosati et al. (2000) from a fit 
to the RDCS data, together with its 2tr errror range. At 
z ~ 0.9 this model lies between a factor of 100 and a factor 
of 1000 below the WARPS XLF measurement. While our 
data thus do not support the RDCS model, they cannot 
strictly rule it out either because of the large uncertain- 
ties in the model parameters. We note though that at 
L x = 5 - 11 x 10 44 h^ 2 erg s" 1 (0.5 - 2.0 keV) the ob- 
served WARPS XLF for very X-ray luminous clusters at 
0.8 < z < 1 is inconsistent with the RDCS prediction at 
greater than 2a significance. 

7. CONCLUSIONS 

The discovery of C1J1226.9+3332, the most X-ray lumi- 
nous distant cluster currently known, adds to the growing 
evidence in favour of an early period of cluster formation 
at redshifts z ^ 1 with little evolution in the cluster abun- 
dance ever since. 

While the small size of the WARPS subsample discussed 
here does not allow us to firmly rule out negative or pos- 
itive evolution at z <~ 1, we stress that our data certainly 
do not require a change in the abundance of very X-ray 
luminous clusters out to the highest redshifts probed by 
current surveys. Specifically, we do not find a decrease in 
the comoving cluster abundance by a factor of more than 
100 (at z <~ 0.9 compared to the local value) as predicted 
by the best-fitting XLF model derived from RDCS data 
(Rosati et al., 2000). 

The lack of significant evolution observed by us at the 
highest redshifts and highest X-ray luminosities is consis- 
tent with the EMSS results at z > 0.8 if the known incom- 
pleteness of the EMSS at these redshifts is corrected for. 

observations in the optical passband for cluster identifications 
0.8 will be presented elsewhere. 
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Our findings also agree with preliminary results from the 
MAssive Cluster Survey (MACS, Ebeling, Edge k. Henry, 
2000) which, boasting greatly improved statistics, finds 
only very mild negative evolution in a measurement of the 
XLF of similarly X-ray luminous clusters at z ~ 0.4. 

Unless cluster evolution is a non-monotonic function of 
redshift and/or X-ray luminosity, our results are in con- 
flict with earlier claims of highly significant (> 3<r) nega- 
tive evolution at lower redshifts (0.3 < z < 0.6) and lower 
luminosities (L x > 3 x 10 44 h^ 2 erg s" 1 , 0.5 - 2.0 keV) 
based on the cluster samples compiled in the EMSS (Gioia 
et al, 1990) and the CfA 160 deg survey (Vikhlinin et al., 
1998a). 



Larger, well defined samples of massive clusters at inter- 
mediate and high redshift as well as spatially resolved clus- 
ter temperatures are needed to actually measure the rate 
of evolution and to place meaningful constraints on the 
cosmological parameters governing structure formation. 
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Table 1 

Positions and heliocentric redshifts of confirmed cluster members 



galaxy R.A. (J2000) Dec (J2000) redshift z 



A 


12 h 


26™ 


58.2 s 


+33° 32' 49" 


0.8911 ±0.0001 


B 


I2 h 


26 m 


57.5 s 


+33° 33' 14" 


0.8794 ± 0.0004 


C 


I2 h 


26 m 


56.1 s 


+33° 32' 23" 


0.8985 ± 0.0005 


D 


12 h 


26 m 


55.6 s 


+33° 32' 13" 


0.9006 ± 0.0005 


E 


12 h 


26 m 


53.0 s 


+33° 32' 50" 


0.8798 ± 0.0003 


F 


12 h 


26 m 


50.6 s 


+33° 32' 46" 


0.8769 ± 0.0004 



Note. — The quoted positions are accurate to better than 
1 arcsec. All redshifts were derived from spectra taken with 
LRIS on Keck-II. The quoted redshift errors are la standard 
deviations from the quoted mean and are based on the ob- 
served wavelengths of individual spectral features. 
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Right Ascension (2000) 



Fig. 1. — I band image of C1J1226. 9+3332 obtained in a four minute exposure with LRIS on Keck-II in February 1999. Overlaid are 
logarithmically spaced contours of the adaptively smoothed X-ray emission as seen with the ROSAT PSPC in the 0.5-2.0 keV band (solid 
lines), and as seen with the ROSAT HRI (dashed lines). The lowest PSPC (HRI) contour lies a factor of 2.5 (1.2) above the background value 
of 2.5 X 10~ 4 (5.4 X 10~ 3 ) ct s" 1 arcmin- 2 ; the levels of adjacent contours differ by a factor of 1.5 (1.2). The FWHM of the PSPC (HRI) 
point-spread function at this off-axis angle (14 arcmin) is 37 (6) arcsec. 
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Fig. 2. — The cumulative X-ray luminosity function (XLF) of very X-ray luminous, distant (z > 0.8) clusters based on C1J0152.7— 1357 and 
C1J1226. 9+3332 as detected in the WARPS survey (solid lines). The strongly shaded region shows the la Poisson error range around the 
XLF assuming a maximal detection rcdshift of z = 1; the hatched region just below delineates the la Poisson error range for the WARPS 
XLF assuming a maximal detection redshift of 2 = 1.4. The local XLF based on the ROSAT Brightest Cluster Sample (Ebeling et al., 1997) 
is shown by the bold, solid line, whereas the prediction of the RDCS XLF evolution model (Rosati et al. 2000) with its 2a error range is 
shown by the dashed line and the lightly shaded region. 



